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We investigate both the intermediate late-time tail and the asymptotic tail behavior of the charged 
\^ , massive Dirac fields in the background of the Kerr-Newman black hole. We find that the intermediate 

' late-time behavior of charged massive Dirac fields is dominated by an inverse power-law decaying 

' tail without any oscillation, which is different from the oscillatory decaying tails of the scalar field. 

. We note that the dumping exponent depends not only on the angular quantum numbers m, the 

^ ' separation constant A and the rotating parameter a, but also on the product seQ of the spin weight 

Q ' of the Dirac field and the charges of the black hole and the fields. We also find that the decay rate 

of the asymptotically late-time tail is t~^^^ , and the oscillation of the tail has the period of 2-k/ ^ 
which is modulated by two types of long-term phase shifts. 

PACS numbers: 04.30.-w, 04.62. -fv, 97.60.Lf. 
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The evolution of field perturbation around a black hole consists roughly of three stages . The first one 
\ is an initial wave burst coming directly from the source of perturbation and is dependent on the initial form 

of the original field perturbation. The second one involves the damped oscillations called the quasinormal 
modes, the frequencies and darnping times of which are entirely fixed by the structure of the background. The 
last one is a power-law tail [d] behavior of the waves at very late time, which is caused by backscattering 
of the gravitational field. 

The late-time evolution of various field perturbations outside a black hole has important irnplications for 
two major aspects of black hole physics: the no-hair theorem and the mass-infiation scenarioQI^]- Therefore, 
since Wheeler 6] 7] introduced the no-hair theorem which states that black holes are completely characterized 
only by three externally observable parameters: mass, electrical charge, and angular momentum, the decay 
' rate of the various fields have been extensively studied [Ij-HJl- Price and Burko 0] studied the massless 

neutral external perturbations and found that the late-time behavior for a fixed r is dominated by the 
factor t~(^'+'^) for each multiple moment. In Refs. the massless late-time tail for the gravitational, 

electromagnetic, neutrino and scalar perturbations had also been considered in the case of the Kerr black 
hole. Starobinskii, Novikov ^| and Burko ^3 analyzed the evolution of a massive scalar field in the Reissner- 
Norstrom background, and they found that, because of the mass term, there are poles in the complex plane 
closer to the real axis than in the massless case, which leads to inverse power-law behavior with smaller indices 
than the massless case. Hod and Piran pointed out that, if the field mass /i is small, namely /lAf ^ 1, 
the oscillatory inverse power-law behavior $ ~ ^-(^+3/2) gin(/it) dominates as the intermediate late-time tails 
in the Reissner-Nordstrom background. Xue and Wang 20] studied the massive scalar wave propagation 
in the background of a Reissner-Nordstrm black hole by using numerical simulations and found that the 
relaxation process depends only on the field parameter and does not depend on the spacetime parameter for 
/iM ^ 1. The very late-time tails of the massive scalar fields in the Schwarzschild and Reissner-Nordstrom 
background were studied in Refs. |21, 22] and in the Kerr case was first investigated numerically in Ref. [2^ . 
and it has been pointed out that the oscillatory inverse power-law behavior of the dominant asymptotic tail 
is approximately given by t"^/^ sin(/zf) which is slower than the intermediate ones. Jing [2J] studied the 
late-time tail behavior of massive Dirac fields in the Schwarzschild black-hole geometry and found that the 
asymptotic behavior of the massive Dirac fields is dominated by a decaying tail without any oscillation in 
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which the dumping exponent depends not only on the multiple number of the wave mode but also on the mass 
of the Dirac field, and also noted that the decay of the massive Dirac field is slower than that of the massive 
scalar field. Recently, late-time tails of self-interacting (massive) scalar fields in a dilaton spacetime and in 
higher dimensional black hole were investigated by Moderski and Rogatko '25*1 "2^ , and they found that the 
intermediate asymptotic behavior of the self-interacting scalar field is determined by an oscillatory inverse 
power-law decaying tail. Konoplya and Molina studied the late-time behavior of massive vector field in 
the background of Schwarzschild black holes, and found three functions with different decay law depending on 
the multiple number I at intermediately late times and a asymptotically late-time decay law cx t~^^^sin{^t) 
which was independent of I. 

Although much attention has been paid to the study of the late-time behaviors of the neutral scalar, 
gravitational, electromagnetic in static and stationary black-hole backgrounds, however, to my best knowledge, 
at the moment only the late-time evolution of the charged massless scalar field was investigated by Hod and 
Piran [25- [sj and their conclusion was that a charged scalar hair outside a charged black hole is dominated 
by a ^-(2^+2) tail which decays slower than a neutral one. Jing |3 recently studied the late-time evolution 
of charged massive Dirac fields in the background of a Reissner-Norstrom black hole note that the dumping 
exponent also depends on the product of the spin weight of the Dirac field and the charges of the black hole 
and the Dirac fields. The main purpose of this paper is extend the study of the late-time tail evolution of the 
charged massive Dirac fields to the stationary Kerr-Newman black-hole background. 

The organization of this paper is as follows. In Sec. II the decoupled charged massive Dirac equations in the 
Kerr-Newman spacetime are presented by using Newman-Penrose formalism. In Sec. Ill the black-hole Green's 
function is introduced and the late-time evolution of the charged massive Dirac fields in the Kerr-Newman 
background is investigated. Sec. IV is devoted to a summary. In appendix we calculate the separation constant 
of the charged massive Dirac equation. 



II. CHARGED MASSIVE DIRAC EQUATIONS IN THE KERR-NEWMAN SPACETIME 

In a curve spacetime the Dirac equations coupled to a electromagnetic field can be expressed as |32| 



V2{V^s-ieA^g)Q^ + i^iPg^0, (2.1) 

where represents the covariant differentiation, and are the two-component spinors representing 
the wave function, Pg is the complex conjugate of Pb, A^^ is a four- vector describes the electromagnetic 
field, and fi and e are the mass and charge of the Dirac particle. Plank units are used, so h = c = G = 1. 
Using the Newman-Penrose formalism ^335 with a null tetrad (1,^,111,171), the equation (|2.1(l becomes 
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For the Kerr-Newman spacetime, the null tetrad can be taken as 

1 -) 



^ ^{{r^ + a"), - A, 0, a), 



= -^{iasinO, 0, 1, -), (2.4) 

with 

A = + + - 2Mr, p^r + iacose, (2.5) 

where M, Q and a represent the mass, charge and angular momentum per unit mass of the Kerr-Newman 
black hole respectively. 
If we take 

po = ^e-*"*+™"^i?_i(r)5_i(0), 
pi ^ e-^"*+""'^i?+i(r)5+i (6*), 
gi ^ e-'^'+'^-^i^+i (r)S'_ 1 (61), 

g° = -ie-*"*+™'*i?_i (r)S'+i(0), (2.6) 

where lo and m are the energy and angular momentum of the Dirac particle, after tedious calculation, the 
equations 1)2.2(1 can be simplified as 

T>nR_ 1 iS'_ 1 H 1 R, 1 S,i ^ —;=iu(r — iacos6)R, i S i , 

2 2 2 +2 T"2 ^2 2 



r — iacos9)R_ i S, i , 

2 '2 ^ 2 2 2 2 2 ~^ 2 

I'o-R iS*,! =£^ii?iiS'_i = —=iti(r + iacos9)R.i S,i , 

2+2 /2 2+2 2 /2 ' + 2 + 2 

iacos9)R_iS_i, (2.7) 
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a iif r - M 

or A A 

+ d iK r-M 

or A A 
9 

-Crt = - -f^ + ncot9, 
o9 

d 

C)n ='^+H + ncot9, 
o9 

K = {r" + a'')LL! — ma — eQr, 

TTl 

H^aujsin9 -. (2.8) 

sin9 

The equation (|2.7|) can be reduced to the following radial and angular parts 

'DoR-i = (A + i^r)i?+i, (2.9) 

AV'' iR^i ^ {\- ip.r)R_i, (2.10) 

-Ci 5+1 = (-A + fiacon9)S_i, (2.11) 

/:^5'_l = {X + ^mcos9)S^i, (2.12) 
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where A is a separation constant. Then we can eliminate S^i (or S'^i) from Eqs. H2.11|l and 12.12|l . and 
obtain the angular equation 



afisinO dSg 
sinO dO \ dO J — 2sA + a^cosO dO 
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(2.13) 



(2.14) 



2sin'^6 2(— 2sA + aficost 

For the slow rotating black hole A^ can be expressed as (see Appendix A for detail) 

A^ ^ \l + Xla + \la^ + ..., 

where Ag, Af and A^ are described by Eqs. IjAQp . (jA10|) and HA14p . 

We can also eliminate R-1/2 (or R+1/2) from Eqs. H2.9|l and H2.10|l . and obtain a radial decoupled equation 
for R+1/2 (or i?_i/2). We find both the decoupled equations for R+1/2 and R-1/2 can be casted into a single 
equation 



with 



^-.iiC^„,^U-^fi + P.R, = 0, (, = ±1) 
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Introducing an usual tortoise coordinate 



dr^, 
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and resolving the equation H2.15|l in the form 
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We will use Eq. ()2.19|l to study the late-time evolution of the charged massive Dirac field. 
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III. LATE-TIME TAIL OF THE CHARGED MASSIVE DIRAC FIELD 



The time evolution of a charged massive Dirac field in Kerr-Newman spacetime is given by 
^s{r.,t) = / [Gs{r,A-,t)dt-^s{<,Q)+dtGs{r,y,-,t)^s{<MK. 



where the black-hole (retarded) Green's function Gs(r*,r^;t) is defined by 



32 32 dH^ 
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with g = ^/r'^ + a? . The causahty condition gives us the initial condition Gs{ri,,r'^]t) = for t <Q. In order 
to get Gs{ri,,r'^;t) we calculate Gs{ri,,r'^;t) through the Fourier transform 



1 



/OQ~\-XC 
Gs{r.,ri;cu)e-''^'duj, (3.3) 
-oo+ic 



where c is a positive constant. The corresponding inversion formula is given by 

/>oo 

Gs{r,,r',;uj) = G,{r,,r',;t)e'^'dt. (3.4) 



The Fourier transform is analytic in the upper half cj-plane and it satisfies the equation 



+ p^^^sj Gs{r,,r:;cj) = S{r, - r'^). (3.5) 



We define auxiliary functions \I'-|-i(rj,, w) and 5'±2(''*,^) which are (linearly independent) solutions to Eq. 
(|2.19|) . Using the solutions 'i'±i and ^±2, the black-hole Greens's function can be constructed as 

G±(r.,r,,^) - I i,^,ir:,c.)^Mr.,c.), > ^^'^^ 

where W±{uj) = VF±(*-t]^, ^±2) = ^±^±2,2: - *±2*±i,x is the Wronskian. 

Morse and Feshbach ^] shown that massive tails exist even in a flat spacetime due to the fact that different 
frequencies forming a massive wave packet have different phase velocities. We will see that at intermediate 
times the backscattering from asymptotically far regions is negligible compared to the flat spacetime massive 
tails that appear here. Hod, Piran, and Leaver fl9, 30] argued that the asymptotic massive tail is associated 
with the existence of a branch cut (in ^2) placed along the interval —fi < uj < ^. This tail arises from the 
integral of the Green function (^(r*, r^; w) around the branch (denoted by G'^(r*, r^; w)) which gives rise to 
an inverse power-law behavior of the Dirac fields. Therefore our goal is to carry out G'^(r*, r!|,; w). 

When both the observer and the initial data are situated far away from the black hole, (a large r or 
equivalently, a low w), we expand the wave equation (|2.19() for the charged massive Dirac field as a power 
series in M/r and Q/r, neglecting terms of order 0(((jj/r)^) and higher, we obtain 

2 2 4Afa)2 _ 2Af^2 _ 2eQw {\/ + 2isM + 2a'm + MIeQ)uj 

— — +L0 - /i -\ 

dr^ r 



A2 _ 2iseQ - Q2(e2 + ^,^^) + (4M2 - a^)^^ 



? = 0, (3.7) 



where ^ = J -p^r^'^ s- We take the transformation 



z = 2\/Jj? — uj^r — 2wr, 

= ^ + \^ -2iseQ-Q^{e^ + ^?) + {AM'^~a^)^x^ 
+ (- + 22sAf + 2am + 8,MeQ)uj, 

k = - - 2Mvo^ (3.8) 



then, Eqs. H3.7(l becomes the confluent hypergeometric equation 

.g + (l + 25-z)f -(i+5-.)<I>^0. (3.9) 
The two basic solutions required in order to build the Green's function can be expressed as 

= v4e-"''(2n7r)i/2+^A/(i + 6 - /c, 1 + 26, 2wr), (3.10) 

*2 = Se-^''(2ti7r)i/2+fc[^( 1 + 6 - fc, 1 + 25, 2wr), (3.11) 
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where A and B are normalization constants. The functions M(k, b, z) and U (fc, 6, z) represent the two standard 
solutions to the confluent hypergeometric equation [s^. U(k^ h, z) is a many- valued function, so there is a cut 
in \E'2. Then, we find that the branch cut which contribute to the Green's function can be expressed as 



G^(r.,r:;i) = 



1 

2^ 



- / F{w)e-''^*duj. 
27r J-n 



W{uj) 



(3.12) 



We can get the following relations with the use of the Eqs. (13.1.32), (13.1.33), and (13.1.34) of Ref. [SE 



^i{2mr) = Ae""^''(2n7r)^+''M(i + 6-fc,l + 26,2n7r), 



+B 



{2-uDrY+^M{- + 6 - fc, 1 + 26, 2vDr) 



^^"^^^ — -e-'^''{2mrY^-Hl(- - 6 - fc, 1 - 26, 2wr). 



Using Eq. (13.1.22) of Ref. we get 



r(26) 



4bm. 



Then, we find that F{'cu) is given by 



26 



M(i + 6 - fc, 1 + 26, 2-cur'-,)M{^ - 6 - fc, 1 - 26, 2cc7r*) 



M(i + 6 + fc, 1 + 26, 2rori)M(i - 6 + fc, 1 - 26, 2ti7r* 



+ 



r(-26)r(i + 6-fc) 
r(26)r(i -6-fc) 



(4^ V^rl) 5+^ 



46cc7 



M(i + 6 - fc, 1 + 26, 2ujri)M{^ + 6 - fc, 1 + 26, 2ror* 



+ e'-^+^''>'''M{^ + 6 + fc, 1 + 26, 2zuri)M{- + 6 + fc, 1 + 26, 2Ti7r* 
In what follows, we will evaluate G"^(r*,r^;t) at the intermediate and very late times, respectively. 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



A. Intermediate behavior of the charged massive Dirac field 

We First focus our attention on the intermediate behavior of the charged massive Dirac field. When the 
late-time behaviors of G'^' at the time scale M ^ r ^ i ^ (MJTj^' which means fc <C 1, following Hod and 

Prian ^3 the effective contribution to the integral in Eq. (|3.3I) is limited to the range vj — 0{y/iJ,/t). 
This is due to the rapidly oscillating term e^*"*, which leads to a mutual cancellation between the positive 
and the negative parts of the integrand. The condition fc <C 1 shows that the 1/r term which describes the 
effect of backscattering off the spacetime curvature from the asymptotically far regions is negligible. 
For the case fc <C 1, noting that M(fc, 6, z) ~ 1 as z ^ 0, we get 

^ ' Amh r(26)r(i - 6) ^ ' ^ 

_ ^ 1 + e(i+2'')" w^^ 

sm(7r6) 21+2&62 r(6)2 ' ' ^^'^"^ 

After tedious calculation, we get 

1 /-M 11 I „(l+26)i7r / r'\h+^ 



4 J_ sm{'Kb) 2^^b^ r(6)2 
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FIG. 1: Graphs of ln\G^ {r^,,ri■,t)\ versus t for different eQ with = 0.01, a = 0.05 and A > 0. For seQ > 0, the 
larger the parameter Q is, the more slowly the perturbation decays. But for seQ < 0, the larger the parameter Q is, 
the more quickly the perturbation decays. 



We work out the integral numerically and present the results in the figures 1 1141 Figure ^ describes 
ln|G'^(r^, r^; t)| versus t for different eQ, and A > 0. We find that the dumping exponent depends on 
the product seQ of the spin weight, the charge of the Dirac particles and the charge of the black hole, and 
seQ > slows the perturbation decay down but seQ < speeds it up. It is the same as A < case. Figure|21ISl 
illustrate ln|G"-^(r^, r'_^;t)\ versus t for different a with Q — 0.8 and ^ = 0.01 and Figure^gives ln\G'~' {r^,,r'^■, t)\ 
versus t for different a with Q = —0.8 and /i = 0.01. Figs. |21 01 and 01 tell us that the dumping exponent 
depends on the quantum number m, the separation constant A and the rotating parameter a, and also show 
that, for both positive and negative electric charge of the black hole, the rotating parameter a slows the decay 
rate down for mX > but speeds it up for mX < 0. 

To compare with the oscillating late-time behavior of the massive scalar field in the Kerr-Newman back- 
ground [3g, we draw the graphs of lii\G^ {ri,,r'^;t)\ versus t with = 0.01, Q — 0.8, a — 0.02 in Fig|Sl It is 
shown that the intermediate late-time behavior of charged massive Dirac field is dominated by a decaying tail 
without any oscillation, and the decay of the massive Dirac field which is affected by the black-hole and field 
parameters is slower than that of the massive scalar field. 

B. Asymptotically late-time tail 

The intermediate tail is not a final pattern of the perturbation decay. At very late times, /it 3> l/{fiM)'^, 
we find 

-fcc.^^^— ^»1, (3.19) 



8 





FIG. 2; Graphs of ln|G'^(r*, rj,; f)| versus t for j — 3/2, with Q = 0.8, fi = 0.01. The rotating parameter a speeds the 
decay up for mX < but slows it down for mX > 0. 




FIG. 3: Graphs of ln\G^ {r^,,ri-t)\ versus t for j = 1/2 and j = 5/2, with Q = 0.8, = 0.01. The rotating parameter 
a speeds the decay up for mX < but slows it down for mX > 0. 



which means that the 1 /r term gives the dominant contribution to the decay. That is to say, the backscattering 
from the asymptotically far regions is important. Using the Eq. (13.5.13) of Ref. we obtain 



F{u;) 



r(l + 2b)r(l-2b)r^< 
26 



r(l + 2b)^T{-2b)Ti} + b - a)ry^ 
26r(26)r(i-6-a) 



J2hi\/ aK )J-2b{ y/aVi, ) - I^2biVar^)I-2b{Vaj\) 
J2b{\/(^)J-2b{^/ar^) + I-2b{y/ari,)I-2b{y/ari,) 



(3.20) 



where a — 8A//i^ — 8eQ^, and I±2b are modified Bessel functions. 

It is easily to know that the factor r(l + 2&)r(l — 2b) /2b in the first term is not a constant. Therefore, we 
use numerical method to study late-time behavior of the first term and present the result in Fig. I^. The 
figure shows, for different value of a, that asymptotically late-time tail arising from the first term is still ~ 
although the factor r(l + 2&)r(l — 2b) /2b is not a constant. 
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FIG. 4: Graphs of ln|G'^(r*, rl,; f)| versus t, with Q = —0.8, /i = 0.01, j = 3/2. The rotating parameter a speeds the 
decay up for m\ < but slows it down for m\ > 0. 




FIG. 5: Graphs of ln|G''^(r*, ; t)\ versus t with fi = 0.01, Q = 0.8, a — 0.02. The dashed lines represent the result of 
the Green function of the massive Dirac field with different j. For comparing, we also show the corresponding result 
of the scalar field with real lines. It is shown that the late-time behavior of massive Dirac fields is dominated by a 
decaying tail without any oscillation, and the decay of the massive Dirac field is slower than that of the massive scalar 
field. 
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FIG. 6: Graphs of ln\G^ {r-,,ri;t)\ versus t with M = 100, = 0.01, Q = O.IM, e = 0.002, s = -1/2 and j = 1/2 
(the left panel is for a = and right one for a — 0.2) for the first term (solid line). The dashed line is ~ log | j|. From 
the figure we know, for different a, that the late-time tail arising from the first term is ~ i. 



Then we consider the contribution of the second term to the decay rate. Because 

r(i + 2&)r(i - 2b) _ r(i + 2b'fT{-2b) 

2& ~ 2bT{2b) ' 

and the factor almost does not affect the asymptotically behavior of the firs term, we define 
r(l + 2b)^r(-2b)r.< 
2br{2b) 

which approximates to a constant. Then, the second term in Eq. H3.20(l can be rewritten as 



A = 



■hbW ar'^)J -2b{\/ ar^,) + I-2b{^/(^)I-2b{^/ar7) 



A 



r(i + 5-fc)^_2.^_, 



2vr r{^-b-k) 



du. 



Using the saddle point integration method, the Eq. (|3.23(l can be written as 



A 
2^ 



where the phase is defined by 



0*0 — 



l + (-l)^''e 



2b„-2i-Kk 



1 + {-lyb^-^i^k ' 

and it remains in the range < </> < 27r, even if k becomes very large. Solving 



— {2nk-ujt) = 0, 
duj 



we get 



Wq 



1/3 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



at very large times (when jit ^ l/(/iAf)^ is satisfied). In this region, most of the waves backscattered by 
the space curvature are cancelled by each other, except for the particular waves with the frequency loq ~ /i. 
Approximating the integral (|3.24(l in the immediate vicinity of ujq, we find 

A^ 



G^(r-*,r:;i) 



2V3 



(27r)5/6(A/A*-eQ)'^'(A*<)"'/' 
. [fit ~ [27r(MM - eQ)]2/3(^t)i/3 „ _ J| 



(3.28) 



The equation shows that the decay rate of the asymptotic late-time tail is exactly t~^^^, and the oscillation has 
the period 27r//i which is modulated by two types of long-term phase shift, the one is [27r(A'//i — e(5)]^/'^(/ii)^/'^ 
represents a monotonously increasing phase shift, the other is 0(wo) represents a periodic phase shift. 

In order to check the above analytical calculation, we also present the numerical result of the second term 
in the Fig. [T] and find that the decay rate of the asymptotically tail is ^ t^^l^ for different value of a. 
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FIG. 7: Graphs of ln|G'^(r^, r^; t)| versus t with M = 100, fi = 0.01, Q = O.IM, e = 0.002, s = -1/2 and j = 1/2 
(the left panel is for a = and right one for a — 0.2) for the second term (solid line). The dashed line is ~ log | | . 
The figure shows, for different a, that the late-time tail arising from the second term is ~ 7576 • 



IV. DISCUSSION AND SUMMARY 



The Teukolsky's master equations for the charged massive Dirac fields in the Kerr-Newman spacetime are 
obtained by means of the Newman-Penrosc formulism. Then, both the intermediate late-time tail and the 
asymptotic tail behavior of the charged massive Dirac fields in the background of the Kerr-Newman black hole 
are studied by considering the intermediate and very late-time solutions of the Teukolsky's master equations. 
The results of the intermediate late-time tail are presented by figures because we can not obtain analytically 
Green's function G'^ (r* , ; t) , as the parameter b in the integrand of the Green's function depends on the 
integral variable ui, especially for a factor \lo/ ^. We learn from the figures that the intermediate late-time 
behavior of charged massive Dirac field is dominated by a decaying tail without any oscillation which is 
different from the oscillatory decaying tails of the scalar field, and the decay of the massive Dirac field is 
slower than that of the massive scalar field. We find that the dumping exponent depends not only on the 
multiple number of the wave mode but also on the product seQ of the spin weight of the Dirac field and the 
charges of the black hole and the fields {seQ < speeds up the decay of the charged massive Dirac fields but 
seQ > slows it down). We also find, for both positive and negative electric charge of the black hole, that 
the rotating parameter a slows the decay rate down for to A > but speeds it up for to A < 0. And at very 
late time, the backscattering off the curvature from the asymptotically far regions is dominant. We find, for 
different value of the angular momentum per unit mass, that the decay rate of the asymptotically late-time 
tail is t^^/^ through both numerical integral method and saddle-point integration method, and the oscillation 
of the tail has the period of 27r//i which is modulated by two types of long-term phase shifts. Our result seems 
to suggest that the oscillatory tail caused by resonance back scattering at asymptotic late time may 

be a quit general feature for the late-time decay of the charged massive Dirac fields in stationary black-hole 
backgrounds. 
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APPENDIX A: THE SEPARATION CONSTANT 



The angular equation H2.13|l can be expressed as 



^2 „d I 



d9^ 



cot9$i 



a^sinO dy 



dO — 2sA + fiacosO dO 



(Al) 
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where 



y = ■\/2sA — aficosO, 
„ nasinO OH 1 nacosO 



A — 2sfiacos6 89 2siv?9 2(— 2sA + a^cosO) 

-H^ + ^cot^e - n^a'^cos^e. (A2) 

For the slow rotating black hole, sCjm and A^ can be expanded as 

A^ = \l + \la + \la^ + 0{a'^), (A3) 

where sPjm are the spherical harmonics with s = ±1/2. The functions s-Pjm and sCjm satisfy the orthogonality 
relations 

sPrm'sP.md^ = I sQ'm'sQmdn = 5,,,6^,^. (A4) 

Using the ordinary perturbation theory, we have 

[Do + Xl]sPjm=0, (A5) 
[^0 + XlhS^ll + [Di + XjUPjm = 0, (A6) 
[Do + Xl]sSl2 + [Di + \l]sS\ll + [D2 + Xl]sPjn^ = 0, (A7) 



with 



(P d 2mscos9 1 1 o 

-C'O = -77^ + COtd— -^r- . + -cot 6 

dO^ dO sm^6 2sm^6 sin^O 4 

/ ficos9 fmi \ 

Di = — — h 2muj — 2sujcos0 , 

4sAo Ao / 



( mfiX^ 3fi^ fiLu 2\ ( -^1 2sm[L\ 
1^-16;^+ Ai^-^j-H^ + ^J'"'' 

AM^-'^^-'-A^osH. (A8) 



IGs^Ao Ao 
From the lowest order equation (|A5|) . we have 

A§ = (J + ^)^ (A9) 
Multiplying equation (|A6|) by sPjm from the left side and integrating it over 0, we obtain 

A? = -2m. + ^ + f - 2..) (AlO) 

Setting 

s'S']m = ^ C^-mS-Pj'm, (All) 



inserting this into equation HA6|) . then multiplying it by gPj'm from the left side and integrating it over 9, we 
have 
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where A = 2s(2'+i) ~ '^slo. It is non-zero only for j' = j ± 1, so we obtain 



2(j + l)2' 



' (j + m + l)(j - m + 1)0- + s + l){j-s + 1) 



(2j + l)(2j + 3) 



1-1 _ _ 



A {j + m) {j -m) {j + s) {j - s) 



(2j-l)(2j + l) 



(A13) 



In order to obtain A2, we multiply equation 1A7|I by gPjm from the left side, and integrate it over 9, then 
we get 



^2 
-^2 



= -2 



J sPjmDlsSjjldil — J sPjmD2sPjmdil 



5^l ^J■uJ 2 2 

los^Ao Ao 



8sAg A2 ; j(j + l) 
1 , 2[3m2-j(j + l)][3.s2-j(j + l)] 



3 3 j(j + l)(2j+3)(2j-l) 



The integrals in the above equations are given by 



' 2 j + 1 , , , , 

'2f~Pl < j,m,l,0|j ,m >< j,-s,l,0|j >, 



2 / 2.7 + 1 



3 V 2j' + 1 



<j,m,2,0\j',m><j,-s,2,0\j',-s > 



(A14) 



(A15) 
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